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Abstract 

The standard Big Bang cosmology predicts the existence of an, as yet undetected, relic 
neutrino background, similar to the photons observed in the cosmic microwave background. 
If neutrinos have mass, then such relic neutrinos are a natural candidate for the dark mat- 
ter of the universe, and indeed were the first particles to be proposed for this role. This 
possibility has however been increasingly constrained by cosmological considerations, partic- 
le | ularly of large-scale structure formation, thus yielding stringent bounds on neutrino masses, 
which have yet to be matched by laboratory experiments. Another probe of relic neutrinos 
l/"") \ is primordial nucleosynthesis which is sensitive to the number of neutrino types (including 
possible sterile species) as well to any lepton asymmetry. Combining such arguments with 
the experimental finding that neutrino mixing angles are large, excludes the possibility of a 
large asymmetry and disfavours new neutrinos beyond those now known. 
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1 Relic neutrinos and the cosmological mass bound 



Several years before neutrinos had even been experimentally detected, Alpher, Follin & Hermann 
[1] noted that they would have been in thermal equilibrium in the early universe ". . . through 
interactions with mesons" at temperatures above 5 MeV; below this temperature the neutrinos 
". . . freeze-in and continue to expand and cool adiabatically as would a pure radiation gas" . 
These authors also observed that the subsequent annihilation of e pairs would heat the photons 
but not the decoupled neutrinos, so by conservation of entropy T u /T would decrease from its 
^ \ high temperature value of unity, down to (4/11) 1 / 3 at T <C m e [1]. 

(Thus the present density of massless relic neutrinos (neglecting a possible chemical potential 
^ . i.e. lepton-antilepton asymmetry) would be 
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where g v = 2 corresponds to left-handed neutrinos and right-handed antineutrinos, and the 
factor 3/4 reflects Fermi versus Bose statistics. This would also be true for massive neutrinos if 
the neutrinos are relativistic at decoupling i.e. for m u <C T^ ec . For a present cosmic microwave 
background (CMB) blackbody temperature To = 2.725 ± 0.002 K [2], the abundance per flavour 
should then be ^ x ^T 3 ~ 111.9 cm . As long as they remain relativistic, the neutrinos 
would retain a Fermi-Dirac distribution with phase-space density 
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since the momentum and temperature would redshift identically.) 

Subsequently, Chiu & Morrison [3] found the rate for e + e~ ^ v e v e in a plasma to be 
T v ~ GpT 5 for the universal Fermi interaction and Zel'dovich [4] equated this to the Hubble 
expansion rate in the radiation-dominated era, 
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H = W — - — , with p = -^9*T , (3) 



where counts the relativistic degrees of freedom, to obtain the 'decoupling' temperature below 
which the neutrinos expand freely without further interactions as 7d ec (^e) — 2MeV. (Neutral 
currents were then unknown so Tf\ cc (v a ) was estimated from the reaction jjl # ez7 e z/„ to be 
12MeV. Later De Graaf [5] noted that these would keep z/^'s coupled to the plasma down 
to the same temperature as zVs. 1 ) However Zel'dovich [4] and Chiu [10] concluded that relic 
neutrinos, although nearly as numerous as the blackbody photons, cannot make an important 
contribution to the cosmological energy density since they are probably massless. 

Interestingly enough, some years earlier Pontecorvo k, Smorodinski [11] had discussed the 
bounds set on the cosmological energy density of MeV energy neutrinos (created e.g. by large- 
scale matter-antimatter annihilation) using data from the Reines-Cowan and Davis experiments. 2 ) 
Not surprisingly these bounds were rather weak so these authors stated somewhat prophetically 
that ". . . it is not possible to exclude a priori the possibility that the neutrino and antineutrino 
energy density in the Universe is comparable to or larger than the average energy density con- 
tained in the proton rest mass" . Zel'dovich & Smorodinski [12] noted that better bounds can be 
set by the limit on the cosmological energy density in any form of matter p m (= 0, m p c ) following 
from the observed present expansion rate Hq and age to of the universe. 3 Of course they were 
still discussing massless neutrinos. Weinberg [13] even speculated whether a degenerate sea of 
relic neutrinos can saturate the cosmological energy density bound and noted that such a sea 
may be detectable by searching for (scattering) events beyond the end-point of the Kurie plot 
in /3-decay experiments. 

Several years later, Gershtem & Zel'dovich [14] made the connection that if relic neutrinos 
are massive, then a bound on the mass follows from simply requiring that 

m u n u <p m . (4) 

Using the general relativistic constraint Q ig Hq < (tt/2) 2 , they derived p m < 2 x 1CU 28 gm cm -3 
(just assuming that ig > 5 Gyr, i.e. that the universe is older than the Earth) and inferred that 
mv e j n^Uft < 400 eV for a present photon temperature of 3 K. Their calculation of the relic 
neutrino abundance was rather approximate — they adopted g u = 4 i.e. assumed massive 
neutrinos to be Dirac particles with fully populated right-handed (RH) states (even though 
they commented en passim that according to the V — A theory such states are non-interacting 
and would thus not be in equilibrium at T& cc ) and moreover they did not allow for the decrease 
in the neutrino temperature relative to photons due to e + e~ annihilation. Nevertheless their 
bound was competitive with the best laboratory bound on m Ue and 10 4 times better than that 
on m v , demonstrating the sensitivity (if not the precision!) of cosmological arguments. 

A better bound of ra V)1 < 130 eV was quoted by Marx & Szalay [15] who numerically 
integrated the cosmological Friedmann equation from decoupling down to the present epoch, 
subject to the condition ig > 4.5 Gyr. Independently Cowsik & McClleland [16] used direct 
limits on f2 m and h to obtain an even more restrictive bound of m v < 8 eV, assuming that 
m v = m Ue = m Vjjl \ however they too assumed that T u = T and that RH states were fully 
populated. As Shapiro, Teukolsky Sz Wasserman [17] first emphasized, even if massive neutrinos 
are Dirac rather than Majorana, the RH states have no gauge interactions so should have 
decoupled much earlier than the left-handed ones. Then subsequent entropy generation by 
massive particle annihilations would have diluted their relic abundance to a negligible level. 4 

1 ln fact Tdec(^M)^T) — 3.5 MeV while Td cc (^e) — 2.3 MeV because of the additional charged current reaction 
[6, 7]. Actually decoupling is not an instantaneous process so the neutrinos are slightly heated by the subsequent 
e + e~ annihilation, increasing the number density (1) by w 1% [8, 9]. 

2 These authors [11] were also the first to suggest searching for high energy neutrinos by looking for upward 
going muons in underground experiments, the basis for today's neutrino telescopes. 

3 The critical density is p c = 3Hq /8tt G n — 1.879 x W~ 29 h 2 gm cm -3 where the Hubble parameter h = 
i/o/100km sec -1 Mpc -1 , so the present Hubble age is Hq 1 = 9.778 ft -1 Gyr. 

4 Although spin-flip scattering (at a rate oc (m„/T) 2 ) can generate RH states, this can be neglected for 
m v <JC 1 MeV. Even if RH neutrinos have new (superweak) interactions, their relic abundance can be no more 
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Thus we arrive at the modern version of the 'Gershtein-Zel'dovich' bound' [20]: the conservative 
limits to > lOGyr and h > 0.4 imply Q m h 2 < 1 i.e. p m < 10.54keVcm~ 3 [21]; combining this 
with the relic neutrino number density, which is ~ 1% larger [9] than in eq.(l), gives: 5 




This is a rather conservative bound since galaxy surveys indicate much tighter constraints on 
the total amount of gravitating (dark) matter in the universe, e.g. the observed 'redshift space 
distortion' [27] suggests Q m « 0.3, averaged over a volume extending several hundred Mpc. 
Together with the Hubble Key Project determination of h = 0.72 ± 0.08 [28], this implies that 
the sum of all neutrino masses cannot exceed about 15 eV. 

Although this has historically been the most restrictive constraint on neutrino masses, it is 
no longer competetive with the direct laboratory bound on the electron neutrino mass from the 
Mainz and Troitsk tritium /3-decay experiments [29]: 

m v < 2.2 eV(95% c.l). (6) 

Although the kinematic mass limits on the other neutrino flavours are much weaker (viz. m v < 
190 keV, m VT < 18.2 MeV [2]), the bound above now applies in fact to all eigenstates [30] given 
the rather small mass-differences indicated by the oscillation interpretation of the Solar (A m? ~ 
7 x 10 -5 eV 2 ) and atmospheric (A m? ~ 3 x 10 -3 eV 2 ) neutrino anomalies [31]. This means that 
the laboratory limit on the sum of neutrino masses (for comparison with the new cosmological 
bounds to be discussed) is presently 6.6 eV; the sensitivity of planned future experiments is at 
the eV level [32]. Moreover relic neutrinos contribute at least Q u h 2 ~ 0.07/93 ~ 8 x 10~ 4 to 
the cosmic budget (assuming a mass hierarchy), nearly as much as visible baryons [33]. 



2 Neutrinos as the galactic 'missing mass' 

The cosmological bound (5) assumes conservatively that neutrinos constitute all of the (dark) 
matter permitted by the dynamics of the universal Hubble expansion. Further constraints must 
be satisfied if they are to cluster on a specified scale (e.g. galactic halos or galaxy clusters) and 
provide the dark matter whose presence is inferred from dynamical measurements. Cowsik & 
McClleland [34] were the first to suggest that neutrinos with a mass of a few eV could naturally 
be the 'missing mass' in clusters of galaxies. This follows from the relation rr? v ~ \ jG\r\M c \ 
(reflecting the Pauli principle) which they obtained by modeling a cluster of mass M c \ as a square 
potential well of core radius r c \ filled with a Fermi-Dirac gas of neutrinos at zero temperature. 
Subsequently Tremaine & Gunn [35] noted that this provides a lower bound on the neutrino 
mass. Although the microscopic phase-space density (2) is conserved for collisionless particles, 
the 'coarse-grained' phase-space density in bound objects can decrease below its maximum value 
of ^/2(27r) 3 during structure formation. Modeling the bound system as an isothermal sphere 
with velocity dispersion a and core radius r\ = 9<7 2 /47r G N p(r c \) then gives 

m„ > 120 eV ( -) ( . (7) 

U00 km sec- 1 / \kpcy v ' 

than ~ 10% of LH neutrinos so the bound on their mass (5) is relaxed to ~ 1 keV [18]. Even if the (post- 
inflationary) universe was not hot enough to bring such interactions into equilibrium, a cosmologically interesting 
abundance can still be generated if the RH states have small mixings with LH states [19] . 

5 If neutrinos are non-relativistic at decoupling, then they drop out of chemical equilibrium with an abundance 
inversely proportional to their self-annihilation cross-section so £l v h 2 ~ (m„/2GeV)~ 2 for m„ <C mz [22, 23]. 
Thus neutrinos with a mass of O(GeV) can also account for the dark matter; however LEP has ruled out 
such (possible 4th generation) neutrinos upto a mass ~ mz/2 [2]. (Conversely Q„h 2 > 1 for the mass range ~ 
100 eV — 2 GeV, which is thus cosmologically forbidden for any stable neutrino having only electroweak interactions 
[24, 25].) For heavier masses upto the highest plausible value of O(TeV), the relic abundance decreases steadily 
due to the increasing annihilation cross-section and remains cosmologically uninteresting [26]. 
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This is indeed consistent with the cosmological upper bound (5) down to the scale of galaxies, 
however there is a conflict for smaller objects, viz. dwarf galaxies which require a minimum 
mass of ~ 100 eV [36, 37]. In fact the central phase space density of observed dark matter cores 
in these structures decreases rapidly with increasing core radius, rather than being constant 
as would be expected for neutrinos [38, 39, 40]. Moreover since neutrinos would cluster more 
efficiently in larger potential wells, there should be a trend of increasing mass-to-light ratio with 
scale. This was indeed claimed to be the case initially [41] but later it was recognised that the 
actual increase is far less than expected [42]. Thus massive neutrinos are now disfavoured as 
the constituent of the 'missing mass' in galaxies and clusters. 

3 The rise and fall of 'hot dark matter' 

Nevertheless such cosmological arguments became of particular interest in the 1980's after the 
ITEP tritium /3-decay experiment claimed a ~ 30 eV mass for the electron neutrino. The 
attention of cosmologists turned to how the large-scale structure (LSS) of galaxies, clusters 
and superclusters would have formed if the universe is indeed dominated by such massive 
neutrinos. The basic picture [43, 44] is that structure grows through gravitational instability 
from primordial density perturbations; these perturbations were first detected by the COsmic 
Background Explorer (COBE) via the temperature fluctuations they induce in the CMB [45]. 
On small scales ( ^ lOMpc) structure formation is complicated by non-linear gravitational 
clustering as well as non-gravitational (gas dynamic) processes but on large scales gravitational 
dynamics is linear and provides a robust probe of the nature of the dark matter. 

Density perturbations in a medium composed of relativistic collisionless particles are subject 
to a form of Landau damping (viz. phase-mixing through free streaming of particles from high 
to low density regions) which effectively erases perturbations on scales smaller than the free- 
streaming length ~ 41 Mpc(m l/ /30eV) _1 [46, 47, 48]. This is essentially the (comoving) distance 
traversed by a neutrino from the Big Bang until it becomes non-relativistic, and corresponds to 
the scale of superclusters of galaxies. Thus huge neutrino condensations (generically in the shape 
of 'pancakes'), containing a mass ~ 3 x lO 15 M (m 1 ,/3OeV)- 2 , would have be gun growing at a 
redshift z cq ~ 7x 10 3 (m^/30eV) when the universe became matter-dominated and gravitational 
instability set in. This is well before the epoch of (re) combination at 

z rec ~ 10^ so the baryons 

were still closely coupled to the photons, while the neutrinos were mildly relativistic (v/c ~ 0.1) 
hence 'hot'. After the universe became neutral, baryonic matter would have accreted into 
these potential wells, forming a thin layer of gas in the central plane of the pancakes. Thus 
superclusters would be the first objects to condense out of the Hubble flow in a 'hot dark matter' 
(HDM) cosmogony, and smaller structures such as galaxies would form only later through the 
fragmentation of the pancakes. 

The gross features of such a 'top-down' model for structure formation are compatible with 
several observed features of LSS, in particular the distinctive 'voids' and 'filaments' seen in large 
galaxy surveys. It was also noted that since primordial density perturbations can begin growing 
earlier than in an purely baryonic universe, their initial amplitude must have been smaller, 
consistent with extant limits on the isotropy of the microwave background. Detailed studies 
[49, 50, 51, 52] found however that galaxies form too late through the breakup of the pancakes, 
at a redshift z <^ 1, counter to observations of galaxies and quasars at z > 4. (Another way of 
saying this is that galaxies should have formed last in an HDM universe, whereas our Galaxy 
is in fact dynamically much older than the local group [53].) There are other difficulties such 
as too high 'peculiar' (i.e. non-Hubble) velocities [54], excessive X-ray emission from baryons 
which accrete onto neutrino clusters [55], and too large voids [56] (although detailed simulations 
[57, 58, 59] showed later that some of these problems had perhaps been exaggerated). 

Therefore cosmologists abandoned HDM and turned, with considerably more success, to cold 
dark matter (CDM) , i.e. particles which were non-relativistic at the epoch of matter-domination 
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[60, 61]. Detailed studies of CDM universes gave excellent agreement with observations of galaxy 
clustering and a 'standard CDM' model for large-scale structure formation was established, 
viz. a critical density CDM dominated universe with an initially scale-invariant spectrum 
of density perturbations [62, 63]. Moreover particle physicists provided plausible candidate 
particles, notably the neutralino in supersymmetric models with conserved .R-parity which 
naturally has a relic abundance of order the critical density [64] . 

4 COBE and the advent of 'mixed dark matter' 

Nevertheless neutrinos were resuscitated some years later as a possible (sub-dominant) compo- 
nent of the dark matter when the CDM cosmogony itself ran into problems. To appreciate the 
background to this it is necessary to recapitulate the essential ingredients of a model for cosmic 
structure formation. A key assumption made concerns the nature of the primordial density 
perturbations which grow through gravitational instability in the dark matter. Cosmologists 
usually assume such fluctuations to have a power spectrum of the scale-free form: 



with n = 1 corresponding to the scale-invariant 'Harrison-Zel'dovich' spectrum. Here 5k = 
j 8p{p_ e -ik-x j g ^ e Fourier transform of spatial fluctuations in the density field (of wavelength 
A = 2ir/k). Moreover the perturbations are assumed to be gaussian (i.e. different phases in the 
plane-wave expansion are uncorrelated) and to be 'adiabatic' (i.e. matter and radiation fluctuate 
together). Concurrent with the above studies concerning the nature of the dark matter, powerful 
support for this conjecture was provided by the development of the 'inflationary universe' model 
[21, 65]. Here the perturbations arise from quantum fluctuations of a scalar field eft, the vacuum 
energy of which drives a period of accelerated expansion in the early universe. The classical 
density perturbation thus generated has a spectrum determined by the 'inflaton' potential V (</>), 
with a power-law index which is dependent on k [66]: 



where M = {SttG^)- 1 / 2 ~ 2.4 x 10 18 GeV is the Planck mass and ★ denotes that this is to 
be evaluated when a mode of wavenumber k crosses the Hubble radius H~ l . For a sufficiently 
'flat potential' (as is necessary to drive enough e-folds of inflation to solve the problems of the 
standard cosmology), the spectrum indeed has n ~ 1, with corrections oc ln(fc). 

Gravitational instability sets in only when the universe becomes matter-dominated and 
this modifies the spectrum on length scales smaller than the Hubble radius at this epoch, 
viz. for k > ~ 80/i _1 Mpc. Thus the characteristics of the dark matter can be encoded 
into a 'transfer function' T{k) which modulates the primordial spectrum; for HDM this is an 
exponentially falling function while for CDM it is a more gradual power-law. Now the power 
spectrum inferred from observations may be compared with theoretical models, but another 
problem arises concerning how we are to normalize the amplitude of the primordial density 
perturbations, 6 particularly since these are in the dark matter and may differ significantly (i.e. 
be 'biased') from the observable fluctuations in the density of visible galaxies. Fortunately, 
the primordial perturbations have another unique observational signature, viz. they induce 
temperature fluctuations in the CMB through the 'Sachs-Wolfe effect' (gravitational red/blue 
shifts) on large angular scales J> 2°, corresponding to spatial scales larger than the Hubble 
radius on the last scattering surface [67]. It was the COBE measurement of these fluctuations 
a decade ago that initiated the modern era of cosmological structure formation studies. 

6 This is determined by the inflaton potential V(4>) but there is, as yet, no 'standard model' of inflation [66]. 



P(k) = (\S k \ 2 )=Ak n , 



(8) 




(9) 
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The quadrupole anisotropy in the CMB measured by COBE [45] allows a determination 
of the fluctuation amplitude at the scale, Hq 1 ~ 3000 h~ l Mpc, corresponding to the present 
Hubble radius. With this normalization it became clear that a Cl v ~ 1 HDM universe indeed 
had too little power on small-scales for adequate galaxy formation. 7 However it also became 
apparent that the 'standard CDM model' when normalized to COBE had too much power on 
small-scales (see Fig. 1). It was thus a logical step to invoke a suitable mixture of CDM and 
HDM to try and match the theoretical power spectrum to the data on galaxy clustering and 
motions [73, 74, 75]. 
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Figure 1: The matter power spectrum as inferred from LSS and CMB data is compared with 
the CDM cosmogony and several variants in this figure reproduced from Scott, White and Silk 
[76]. We see top left that the excess small-scale power in the COBE-normalized 'standard' CDM 
model with scale-invariant adiabatic fluctuations (n = 1,O b = 0.03, h = 0.5), is reduced in the 
MDM model which has a substantial neutrino component (0„ = 0.3). This can also be achieved 
by 'tilting' the spectrum to n = 0.9 as shown top right (f2 B = 0.1, h = 0.45). The middle panels 
show the expectations in open universes, while the bottom panels correspond to a flat universe 
with a cosmological constant (the bottom right figure assumes a significant gravitational wave 
contribution to the COBE signal). 



In fact the possibility that the dark matter may have both a hot and a cold component had 
been discussed several years earlier by Shafi & Stecker, motivated by theoretical considerations 
of SUSY GUTs [77] , and the possibilities for solving the problems of the standard CDM model 
had been noted [78]. In the post-COBE era, a number of detailed studies of mixed dark matter 

7 To save HDM would require new sources of small-scale fluctuations, e.g. relic topological defects [68, 69, 70], 
or isocurvature primordial perturbations [71, 72] — to date however there is no evidence for either. 
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(MDM) universes were performed and a neutrino fraction of about 20% was found to give the 
best match with observations [79, 80, 81, 82, 83, 84]. The implied neutrino mass was ~ 5eV, 
presumably that of the v r given the usual hierarchy implied by the see-saw mechanism [85]. 
The suppression of small-scale power delays the epoch of galaxy formation, so constraints on 
the HDM component can also be obtained from the abundance of high redshift objects such as 
QSOs and Ly-a systems; this gave an upper bound of 4.7 eV on the neutrino mass [86]. More 
baroque schemes in which two neutrinos have comparable masses (m v ~ m Ur ~ 2.5 eV) were 
also constructed [87] seeking to reconcile the LSND report of neutrino oscillations [88] with the 
atmospheric neutrino anomaly which had just been reported by Kamiokande [89]. This was an 
exciting time for neutrino cosmology as both laboratory data and astronomical observations 
supported the possibility that a substantial fraction of the cosmological mass density is in the 
form of massive neutrinos. 

However another way to reconcile a CDM universe with the small-scale observations is to 
relax the underlying assumption that the primordial spectrum is strictly scale-invariant. As 
shown in Fig.l, a 'tilted' spectrum with P{k) oc /c ' 9 also gives a good fit to the data [90]. 
At first sight this might strike one as simply introducing an additional parameter (although 
this is arguably no worse than introducing an additional form of dark matter). However one 
should really ask why the spectrum should be assumed to have a power-law index n = 1 in 
the first place. As indicated in Eq.(9), n in fact varies slowly with k and is determined by the 
slope and curvature of the scalar potential at the epoch when the fluctuation at a specified 
value of k crosses the Hubble radius. The corresponding number of e-folds before the end of 
inflation is just N+(k) ~ 51 + ln ( 3000 fe-i m pc ) ' f° r typical choices of the inflationary scale, reheat 
temperature etc. We see that fluctuations on the scales (~ 1 — 3000 Mpc) probed by LSS and 
CMB observations are generated just 40 — 50 e-folds before the end of inflation. It is quite 
natural to expect the inflaton potential to begin curving significantly as the end of inflation 
is approached, especially in 'new inflation' (small field) models. There are certainly attractive 
models of inflation in which the spectrum is significantly tilted in this region, in particular an 
inflationary model based on N = 1 supergravity [91] naturally gives n(k) ~ (N* — 2)/(N ir + 2) ~ 
0.9 at these scales. With such a scale- dependent tilt for the primordial spctrum, the LSS data 
can be fitted reasonably well, with no need for any HDM component [92]. (Of course in the 
absence of a 'standard model' of inflation, it might be argued that the inflationary spectrum 
may instead have n > 1 thus allowing a larger HDM component.) 

As Fig.l shows, yet another way of suppressing small-scale power in the CDM cosmogony is 
to decrease the matter content of the universe, since this postpones the epoch of matter-radiation 
equality and thus shifts the peak of power spectrum to larger scales. Furthermore the spatial 
geometry can be maintained flat if there is a cosmological constant with = 1 — £l m ~ 0.7. 
Evidence for just such a cosmology (ACDM) has come subsequently from observations of the 
Hubble diagram of Type la supernovae which suggest that the expansion is in fact accelerating 
[93, 94], coupled with the observation that O m does not exceed ~ 0.3 even on the largest 
scales probed [95]. Assuming such a ACDM cosmology and a scale-invariant power spectrum 
gives a limit of f v = fi„/fi m < 0.13 (95% c.l.) using the power spectrum of galaxy clustering 
determined from the Two degree Field galaxy redshift survey (2dFGRS) as shown in Fig. 2. 
This corresponds to an upper bound of J2 m v < 1-8 eV(95% c.l.) on the sum of neutrino masses, 
adopting 'concordance' values of Q m and h [97]. If the spectral index is allowed to vary in the 
range n = 1 ± 0.1, then this bound is relaxed to 2.1 eV. 

5 CMB anisotropy and limits on neutrino dark matter 

It is clear that there are uncertainties in the determination of the HDM component from LSS 
data alone. Fortunately it is possible to reduce these substantially by examining an independent 
probe of the primordial power spectrum, viz. temperature fluctuations in the CMB. 
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Figure 2: The 2dFGRS power spectrum of galaxy clustering is compared with the expectations 
for no neutrino dark matter (solid line), Q u = 0.01 (dashed line) and Q. v = 0.05 (dot-dashed 
line). A scale- invariant spectrum is assumed and the cosmological parameters adopted are 
n m = 0.3, ft A = 0.7, h = 0.7, Q B h 2 = 0.02 (from Elgaroy et al. [97]). 

In general a skymap of the CMB temperature can be decomposed into spherical harmonics 



where the I th multipole corresponds to an angle 8° ~ 200/7 and probes spatial scales around 
k~ l ~ 6000/i -1 / -1 Mpc. In inflationary theories, the fluctuations are gaussian so the co-efficients 
a™ are independent stochastic variables with zero mean and variance C/ = (|a™| 2 ); each C\ has 
a x 2 distribution with (21 + 1) degrees of freedom [60, 98]. For an assumed set of cosmological 
parameters and given the primordial density perturbation spectrum, the C/'s can be determined 
by solution of the Einstein-Boltzmann equations which describe how the different components 
(photons, ions, electrons, neutrinos . . . ) evolve [99, 100, 101]. Thus theoretical estimates of the 
power at each multipole can be compared with observations. The low multipoles (large spatial 
scales) are sensitive to the primordial spectrum alone, 8 but the measurements in the region are 
particularly uncertain, both because of uncertainties in the foreground substraction and also 
because there are fewer independent measurements for low multipoles ('cosmic variance'). For 
example, by measuring the first ~20 multipoles COBE could only fix n = 1.2 ±0.3 [45] so could 
not discriminate between a scale-invariant and a mildly tilted spectrum. However subsequent 
ground-based experiments with angular resolution far superior to COBE's have now measured 
the power at higher multipoles [103]. The dominant features in the power spectrum here are the 
'acoustic peaks', the most prominent being at / ~ 200, arising from oscillations of the coupled 
plasma-photon fluids at last scattering [104]. The position of the first peak is a measure of the 
horizon length at the epoch of (re)combination of the primordial plasma and thus provides a 
measure of space curvature — observations by the BOOMERanG [105] and MAXIMA [106] 
experiments were the first to show that this is in fact close to zero. Taken together with the 
earlier recognition that CDM alone cannot make up the critical density, this has led to the 
widespread adoption of the ACDM model in which fi m ~ 0.3, ~ 0.7 [95, 96]. 

As seen in Fig. 3, the expectations for CMB anisotropy in a MDM universe do not differ 
significantly from a CDM universe having the same initial perturbation spectrum. Thus it 

8 In principle, primordial gravitational waves can also make a contribution here but this is expected to be 
negligible in 'small- field' inflationary models [91, 102]. 
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Figure 3: The angular power spectrum of CMB anisotropy (assuming a scale-invariant spec- 
trum) in CDM and MDM universes, as calculated by Dodelson, Gates and Stebbins [107]. 

is clear that by combining CMB and LSS data, it would be possible to determine whether 
the suppression of small-scale power is intrinsic to the primordial spectrum of inflationary 
fluctuations, or is induced by a HDM component. However there are additional uncertainties 
in the input values of the other cosmological parameters, in particular the values of the Hubble 
parameter and of the baryon density which is usually inferred from considerations of primordial 
nucleosynthesis. In analysing the observational data, various 'priors' are often assumed for 
these quantities on the basis of other observations. For example a detailed likelihood analysis 
[108] yielded the neutrino density fraction f u to be at most 0.08 from CMB observations alone, 
decreasing to 0.06 if LSS data from the PSCz survey is added, and further to 0.04 if the value of 
h = 0.72 ±0.08 [28] is adopted. Another analysis [109] used the more precise 2dFGRS data [97], 
together with additional constraints on the matter density from the SNIa data (O m = 0.28±0.14 
assuming a flat universe) [93, 94], and on the baryon density from primordial nucleosynthesis 
(Q B h 2 = 0.02±0.002 (95% c.l.) [Ill]); using the analytic result for the suppression of the power 
spectrum, AP/P ~ — 8/„ [110], a bound of m v < 2.5 — 3 eV was obtained. 




Figure 4: The marginalized cumulative probability of Q. v h 2 is shown, based on a fit to the 
WMAP data on CMB anisotropies, together with the 2DFGRS data on of galaxy clustering 
(from Spergel et al. [113], courtesey of WMAP Science Team). 
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It had been estimated [110] that the CMB data expected from the Wilkinson Microwave 
Anisotropy Probe (WMAP) [112], combined with LSS data from the Sloane Digital Sky Survey 
(SDSS) [118], will provide sensitivity to neutrino mass at the eV level. Indeed the recently 
announced first results from WMAP, combined with 2dFGRS, have already set the bound 
n u h 2 < 0.0076 corresponding to < 0.7 eV(95% c.L), as shown in Fig. 4 [113]. This 

severely restricts [114], but does not altogether rule out [115], a fourth (sterile) neutrino with a 
mass of 0(1) eV as suggested by the LSND experiment [88]. Combined with the new KamLAND 
data, this also restricts the accessible range for the observation of neutrinoless (3(3 decay [116], 
which has been claimed to have been seen already implying an effective Majorana mass of 
\{m v )\ = 0.39t°; 4 |eV (95% c.L) [117]. 



6 The nucleosynthesis limit on N v 

Hoyle & Taylor [119] as well as Peebles [120] had emphasized many years ago that new types 
of neutrinos (beyond the v e and then known) would boost the relativistic energy density 
hence the expansion rate (3) during Big Bang nucleosynthesis (BBN), thus increasing the yield 
of 4 He. Shvartsman [121] noted that new superweakly interacting particles would have a similar 
effect. Subsequently this argument was refined quantitatively by Steigman, Schramm & Gunn 
[122]. In the pre-LEP era when the laboratory bound on the number of neutrino species was 
not very restrictive [123], the BBN constraint was used to argue that at most one new family 
was allowed [124, 125], albeit with considerable uncertainties [126]. Although LEP now finds 
N v = 2.994 ±0.012 [2], the cosmological bound is still important since it is sensitive to any new 
light particle [127], not just S77(2) L doublet neutrinos, so is a particularly valuable probe of 
new physics, e.g. neutrinos coupled to new gauge bosons expected in string models [126, 128]. 

BBN limits on neutrinos come mainly from the observational bounds on the primordial 
4 He abundance, termed Y p by astronomers. This is proportional to the neutron fraction which 
'freezes out' at n/p ~ exp[— (m n — m p )/Tf r ] when the weak interaction rate (oc G*|T 5 ) falls 
behind the Hubble expansion rate (3), at T fr ~ (g*G N /G^) 1 / 6 « 1 MeV. The presence of 
additional neutrino flavors (or of any other relativistic species) at this time increases hence 
the expansion rate H, leading to a larger value of Tf r , n/p, and thus Y p . In the Standard Model, 
the number of relativistic particle species at 1 MeV is 

^ = 5.5 + 1^, (11) 

where 5.5 accounts for photons and e^, and N v is the number of (massless) neutrino fla- 
vors. (The energy density of new light fermions i is equivalent to an effective number A N v = 
E i (<7i/2)(Tj/T l/ ) 4 of additional doublet neutrinos, where Ti/T u follows from considerations of 
their (earlier) decoupling.) In Fig. 5, we show the expected abundance of 4 He taking N v = 3, 
as a function of the density of baryons normalized to the blackbody photon density: 

r/ = — ~ 2.728 x lO" 8 ^ 2 , (12) 

n 7 

with 7710 = n/W~ 10 . The computed 4 He abundance scales linearly with AN U but it also 
increases logarithmically with n [129, 130]: 



Ay ~ O.O^AiVV + 0.01 In 




Thus to obtain a bound on N u requires an upper limit on Y p and a lower limit on rj. The latter 
is poorly determined from direct observations of luminous matter [33] so must be derived from 
the abundances of the other synthesized light elements, D, 3 He and Li, which are power-law 
functions of n as seen in Fig. 5. The complication is that these abundances are substantially 
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altered in a non-trivial manner during the chemical evolution of the galaxy, unlike Yp( 4 He) 
which just increases by a few percent due to stellar production. (This can be tagged via the 
correlated abundance of oxygen and nitrogen which are made only in stars.) Even so, chemical 
evolution arguments have been used to limit the primordial abundances of D and 3 He and 
thus derive increasingly severe bounds on N v [131, 132, 133], culminating in one below 3 [134]! 
However a more conservative view [135] is that there is no 'crisis' with BBN if we recognize 
that chemical evolution arguments are unreliable and consider only direct measurements of light 
element abundances. 



0.005 



0.01 



Q. B h 2 



0.02 



0.03 
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Baryon-to-photon ratio r\ 1Q 



7 8 9 10 



Figure 5: The primordial abundances of 4 He, D, 3 He and 7 Li as predicted by the standard 
BBN model compared to observations — smaller boxes: 2a statistical errors; larger boxes: ±2cr 
statistical and systematic errors added in quadrature (from Fields and Sarkar [2]). 



Fig. 5 shows the current status of such measurements, as reviewed in more detail elsewhere 
[2, 136]. We observe 4 He in clouds of ionized hydrogen (H II regions), the most metal-poor of 
which are in dwarf blue compact galaxies (BCGs). There is now a large body of data on 4 He and 
C, N, O in these systems which confirm that the small stellar contribution to helium is positively 
correlated with 'metal' production; extrapolating to zero metallicity gives the primordial 4 He 
abundance [137]: 

Y p = 0.238 ± 0.002 ± 0.005 . (14) 
Here and subsequently, the first error is statistical, and the second is an estimate of the system- 
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atic uncertainty. The latter clearly dominates, and is based on the scatter in different analyses 
of the physical properties of the H II regions [138, 139]. Other extrapolations to zero metal- 
licity give Y p = 0.244 ± 0.002 [140], and Y p = 0.235 ± 0.003 [141], while the average in the 5 
most metal-poor objects is Y p = 0.238 ± 0.003 [142]. The value in Eq.(14), shown in Fig. 5, is 
consistent with all these determinations. 

The systems best suited for Li observations are hot, metal-poor stars belonging to the 
halo population (Pop II) of our Galaxy. Observations have long shown that Li does not vary 
significantly in such stars having metallicities less than 1/30 of Solar — the 'Spite plateau'. 
Recent precision data suggest a small but significant correlation between Li and Fe [143] which 
can be understood as the result of Li production from cosmic rays [144]. Extrapolating to zero 
metallicity one arrives at a primordial value [145] 

Li/H| p = (1.23 ± 0.06lg;l +0 - 56 ) x 10~ 10 . (15) 

The last error is our estimate of the maximum upward correction necessary to allow for possible 
destruction of Li in Pop II stars, due to e.g. mixing of the outer layers with the hotter interior 
[146, 147]. Such processes can be constrained by the absence of significant scatter in the Li-Fe 
correlation plot [143], and through observations of the even more fragile isotope 6 Li [144]. 

In recent years, high-resolution spectra have revealed the presence of D in quasar absorption 
systems (QAS) at high-redshift, via its isotope-shifted Lyman-a absorption. It is believed 
that there are no astrophysical sources of deuterium, so any measurement of D/H provides a 
lower limit to the primordial abundance and thus an upper limit on 77; for example, the local 
interstellar value of D / H = (1.5 ± 0.1) x 10~ 5 [148] requires r?i < 9. Early reports of D / H > 
10~ 4 towards 2 quasars (Q0014+813 [149] and PG1718+4807 [150]) have been undermined by 
later analyses [151, 152]. Three high quality observations yield D/H = (3.2 ± 0.3) x 10~ 5 
(PKS1937-1009), (4.0 ±0.7) x 10~ 5 (Q1009+2956), and (2.5±0.2) x 10~ 5 (HS0105+1619); their 
average value 

D/H = (3.0 ±0.4) x 10" 5 (16) 

has been widely promoted as the primordial abundance [153]. Recently, the same group [154] 
have provided another measurement: D/H = (2.42t°;||) x 10~ 5 (HS 243+3057). However the 
observed dispersion in the measurements suggests either that systematic uncertainties have been 
underestimated, or that there is intrinsic dispersion in the D abundance in QAS. Other values 
have been reported in different (damped Lyman-a) systems which have a higher column density 
of neutral H, viz. D / H = (2.24±0.67) x 10~ 5 (Q0347-3819) [155] and D / H = (1.65±0.35) x 10~ 5 
(Q2206-199) [156]. Moreover, allowing for a more complex velocity structure than assumed in 
these analyses raises the inferred abundance by upto ~ 50% [157]. Even the ISM value of D/H 
now shows unexpected scatter of a factor of 2 [158]. All this may indicate significant processing 
of the D abundance even at high redshift. Given these uncertainties, it is prudent to bound 
the primordial abundance with an upper limit set by the non-detection of D absorption in a 
high-redshift system (Q0130-4021) [159], and the lower limit set by the local interstellar value 
[148], both at 2a: 

1.3 x 10~ 5 < D/H| P < 9.7 x 10~ 5 . (17) 

For 3 He, the only observations available are in the Solar system and (high-metallicity) H II 
regions in our Galaxy [160]. This makes inference of the primordial abundance difficult, a 
problem compounded by the fact that stellar nucleosynthesis models for 3 He are in conflict 
with observations [161]. Such conflicts can perhaps be resolved if a large fraction of low mass 
stars destroy 3 He by internal mixing driven by stellar rotation, consistent with the observed 
12 C/ 13 C ratios [162]. The observed abundance 'plateau' in H II regions then implies a limit on 
the primordial value of 3 He/H < (1.9 ± 0.6) x 10~ 5 [160], which is consistent with the other 
abundance constraints we discuss. 

The overlap in the r/ ranges spanned by the larger boxes in Fig. 5 indicates overall concor- 
dance between the various abundance determinations. Accounting for theoretical uncertainties 
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as well as the statistical and systematic errors in observations, there is acceptable agreement 
among the abundances when [2] 

2.6 < 7710 < 6.2 . (18) 

However the agreement is far less satisfactory if we use only the quoted statistical errors in the 
observations. As seen in Fig. 5, 4 He and 7 Li are consistent with each other but favor a value 
of r\ which is discrepant by ~ 2a from the value r/io = 5.9 ± 0.4 indicated by the D abundance 
(16) alone. It is important to note that it is the latter which essentially provides the widely- 
quoted 'precision ' determination of the baryon density: Q B h 2 = 0.02 ± 0.002 (95% c.l.) [111]. 
Additional studies are required to clarify if this discrepancy is real. 

The recent WMAP data on CMB fluctuations implies a baryon abundance of r/io = 6.5+03 if 
the primordial perturbations are assumed to have a power-law form [113]; the value decreases by 
~ 9% if the spectral index is allowed to vary with scale, as is now indicated by the observations. 
Although broadly consistent with the determinations based on primordial abundances, it is 
clear that this exacerbates the tension with the 4 He and 7 Li measurements mentioned above. 
It has been proposed that the CMB determination of 7710 can be used as an input into BBN 
calculations [163]. However it would be advisable to await further LSS data from SDSS [118] to 
pin down more precisely the primordial fluctuation spectrum, to which such determinations are 
sensitive. For example, allowing for a 'step' in the primordial spectrum at k ~ 0.05/i Mpc -1 , as 
was indicated by data from the APM galaxy survey [164, 165], can decrease the baryon density 
inferred from the WMAP+2DFGRS data by upto ~ 40% [166]. 

Given this situation, it is still neccessary to be conservative in evaluating bounds on N v 
from BBN. An analysis based on simple \ 2 statistics and taking into account the correlated 
uncertainties of the elemental yields, gives [167] 

2 < N v < 4 (95% c.l). (19) 

Tighter bounds can of course be obtained under less conservative assumptions, e.g. adopting the 
D abundance (16) requires N v < 3.2 [168]. If true, such bounds rule out the possibility of 'sterile' 
neutrinos since these would be brought into equilibrium through mixing with the active neutri- 
nos in all currently viable schemes [169, 170]. The WMAP determination of the baryon density 
seems to require N v less than 3 [171] but is consistent with the Standard Model taking system- 
atic uncertainties in the elemental abundances into account [114, 172]. Moreover as emphasized 
above, the WMAP determination of r\ will have to ve lowered if the primordial spectrum is not 
scale-free and this will considerably ease the present tension. (The WMAP+2dFGRS data by 
themselves restrict 1.8 < N v < 5.7 (95% c.l.) assuming a flat universe [173, 174].) 

The above limits hold for the standard BBN model and one can ask to what extent they 
can be modified if plausible changes are made to the model. For example, there may be an 
initial excess of electron neutrinos over antineutrinos, parametrised by a 'chemical potential' 
£ e in the distribution function: f u oc [exp(p/T — £ e ) + 1] _1 . Then n — p equilibrium is shifted 
in favour of less neutrons, thus reducing the 4 He abundance by a factor rj exp(— £ e ) [129]. 
However the accompanying increase in the relativistic energy density speeds up the expansion 
rate and increases the n/p ratio at freeze-out, leading to more 4 He, although this effect is 
smaller. For neutrinos of other flavours which do not participate in nuclear reactions, only the 
latter effect was presumed to operate, allowing the possibility of balancing a small chemical 
potential in v e by a much larger chemical potential in v^ T , and thus substantially enlarging 
the concordance range of 7710 [175]. However the recent recognition from Solar and atmospheric 
neutrino experiments that the different flavours are maximally mixed, no longer permits such 
a hierarchy of chemical potentials [176, 177, 178, 179], thus ruling out this possible loophole. 
Consequently the relic neutrino abundance cannot be significantly different from the value (1). 
A small chemical potential in electron neutrinos of £ e 0.08 is however still possible and could 
be another explanation for why the 4 He abundance seems to be lower than the value that would 
be expected on the basis of the measurements of the D abundance or the CMB anisotropy. 
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Another possible change to standard BBN is to allow inhomogeneities in the baryon distri- 
bution, created e.g. during the QCD (de) confinement transition. If the characteristic inhomo- 
geneity scale exceeds the neutron diffusion scale during BBN, then increasing the average value 
of 7] increases the synthesised abundances such that the observational limits essentially rule 
out such inhomogeneities. However fluctuations in 77 on smaller scales will result in neutrons 
escaping from the high density regions leading to spatial variations in the n/p ratio which might 
allow the upper limit to n to be raised substantially [180]. Recent calculations show that D 
and 4 He can indeed be matched even when n is raised by a factor of ~ 2 by suitably tuning the 
amplitude and scale of the fluctuations, but this results in unacceptable overproduction of 7 Li 
[181]. A variant on the above possibility is to allow for regions of antimatter which annihilate 
during or even after BBN; however the 7 Li abundance again restricts the possibility of raising 
the limit on r\ substantially [182]. Finally the synthesised abundances can be altered if a relic 
massive particle decays during or after BBN generating electromagnetic and hadronic showers 
in the radiation-dominated plasma. Interestingly enough the processed yields of D, 4 He and 
7 Li can then be made to match the observations even for a universe closed by baryons [183], 
however the production of 6 Li is excessive and argues against this possibility. 

In summary, standard BBN appears to be reasonably robust (although non-standard pos- 
sibilities cannot be definitively ruled out) and consistent with N v = 3, leaving little room for 
new physics. Systematic uncertainties in the elemental abundance determinations need to be 
substantially reduced in order to make further progress. 

7 Conclusions 

The study of neutrinos in cosmology has had a long history. However as both theory and 
observations have improved, the fascinating possibilities that had been raised in early work 
have gradually been eliminated. In particular neutrino (hot) dark matter is no longer required 
by our present understanding of large-scale structure. Cosmological bounds on neutrino masses 
continue to be more restrictive than laboratory limits and all known species are now required to 
be sufficiently light that the rich phenomenology of unstable neutrinos is now largely of historical 
interest. Moreover there is no motivation from considerations of primordial nucleosynthesis 
for invoking new sterile neutrinos. Nevertheless such neutrinos can exist (e.g. in the 'bulk' 
[184, 185]) and there may yet be surprises from the Big Bang in store for us. 

References 

[1] R.A. Alpher, J.W. Follin & R.C. Hermann, Phys. Rev. 92 (1953) 1347. 

[2] K. Hagiwara et al. [Particle Data Group Collaboration], Phys. Rev. D 66 (2002) 010001. 

[3] H.Y. Chiu & P. Morrison, Phys. Rev. Lett. 5 (1960) 573. 

[4] Ya.B. Zel'dovich, Adv. Astron. Astrophys. 3 (1965) 241; Sov. Phys. Usp. 9 (1967) 602. 
[5] T. De Graaf, Lett. Nuovo Cim. 4 (1970) 638. 

[6] D.A. Dicus, E.W. Kolb, A.M. Gleeson, E.C. Sudarshan, V.L. Teplitz & M.S. Turner, Phys. 
Rev. D 26 (1982) 2694. 

[7] K. Enqvist, K. Kainulainen & V. Semikoz, Nucl. Phys. B 374 (1992) 392. 

[8] S. Dodelson & M.S. Turner, Phys. Rev. D 46 (1992) 3372. 

[9] S. Hannestad & J. Madsen, Phys. Rev. D 52 (1995) 1764. 



14 



[10] H. Chiu, Ann. Rev. Nucl. Sci. 16 (1966) 591. 

[11] B. Pontecorvo k Ya. Smorodinski, Sov. Phys. JETP 14 (1962) 173. 

[12] Ya.B. Zcl'dovich k Ya. Smorodinski, Sov. Phys. JETP 14 (1962) 647. 

[13] S. Weinberg, Phys. Rev. 128 (1962) 1457. 

[14] S.S. Gershtein k Ya.B. Zeldovich, JETP Lett. 4 (1966) 120. 

[15] G. Marx k A.S. Szalay, Proc. Neutrino'72, Balatonfiird, Vol.1 (1972) p.123; 
A.S. Szalay k G. Marx, Astron. Astrophys. 49 (1976) 437. 

[16] R. Cowsik k J. McClelland, Phys. Rev. Lett. 29 (1972) 669. 

[17] S.L. Shapiro, S.A. Teukolsky k I. Wasserman, Phys. Rev. Lett. 45 (1980) 669. 

[18] K.A. Olive k M.S. Turner, Phys. Rev. D 25 (1982) 213. 

[19] A.D. Dolgov k S.H. Hansen, Astropart. Phys. 16 (2002) 339. 

[20] J. Bernstein k G. Feinberg, Phys. Lett. B 101 (1981) 39 [Erratmn-ibid. B 103 (1981) 470]. 

[21] E.W. Kolb k M.S. Turner, The Early Universe (Addison- Wesley, 1990). 

[22] B.W. Lee k S. Weinberg, Phys. Rev. Lett. 39 (1977) 165. 

[23] M.I. Vysotsky, A.D. Dolgov k Ya.B. Zeldovich, JETP Lett. 26 (1977) 188. 

[24] K. Sato k M. Kobayashi, Prog. Theor. Phys. 58 (1977) 1775. 

[25] D.A. Dicus, E.W. Kolb k V.L. Teplitz, Astrophys. J. 221 (1978) 327. 

[26] K. Enqvist, K. Kainulainen k J. Maalampi, Nucl. Phys. B 317 (1989) 647. 

[27] E. Hawkins et al, arXiv:astro-ph/0212375. 

[28] W.L. Freedman et al. [Hubble Key Project Collaboration], Astrophys. J. 553 (2001) 47. 
[29] C. Weinheimer, arXiv:hep-ex/0210050. 

[30] V.D. Barger, T.J. Weiler k K. Whisnant, Phys. Lett. B 442 (1998) 255. 

[31] M.C. Gonzalez-Garcia k Y. Mr, arXiv:hep-ph/0202058. 

[32] A. Osipowicz et al. [KATRIN Collaboration], arXiv:hep-ex/0109033. 

[33] M. Fukugita, C.J. Hogan k P.J.E. Peebles, Astrophys. J. 503 (1998) 518 

[34] R. Cowsik k J. McClleland, Astrophys. J. 180 (1972) 7. 

[35] S. Tremaine k J.E. Gunn, Phys. Rev. Lett. 42 (1979) 407. 

[36] D.N. Spergel, D.H. Weinberg k J.R. Gott, Phys. Rev. D 38 (1988) 2014. 

[37] G. Lake, Astron. J. 98 (1989) 1253. 

[38] A. Burkert, in Aspects of Dark Matter in Astro and Particle Physics, eds H.V. Klapdor k 
Y. Ramachers (World Scientific, 1997) p.35. 

[39] P. Salucci k A. Sinibaldi, Astron. Astrophys. 323 (1997) 1. 

15 



[40] J.A. Sellwood, Astrophys. J. 540 (2000) LI. 

[41] D.N. Schramm k G. Steigman, Astrophys. J. 243 (1981) 1. 

[42] G.R. Blumenthal, S.M. Faber, J.R. Primack & M.J. Rees, Nature 311 (1984) 517. 
[43] T. Padmanabhan, Structure Formation in the Universe (Cambridge Univ Press, 1993). 
[44] P.J.E. Peebles, Principles Of Physical Cosmology (Princeton Univ Press, 1993). 
[45] G.F. Smoot et al, Astrophys. J. 396 (1992) LI. 

[46] A.G. Doroshkevich, Ya.B. Zeldovich, R.A. Sunyaev & M.Yu. Khlopov, Sov. Astron. Lett. 
6 (1980) 251, 265 

[47] J.R. Bond, G. Efstathiou & J. Silk, Phys. Rev. Lett. 45 (1980) 1980. 

[48] A.G. Doroshkevich, M.Yu. Khlopov, R.A. Sunyaev, Ya.B. Zeldovich & A.S. Szalay, Ann. 
N.Y. Acad. Sci. 375 (1981) 32. 

[49] P.J.E. Peebles, Astrophys. J. 258 (1982) 415. 

[50] A. A. Klypin & S.F. Shandarin, Mon. Not. R. Astr. Soc. 204 (1983) 891. 
[51] S.D.M. White, C.S. Frenk & M. Davis, Astrophys. J. 274 (1983) LI. 
[52] J.R. Bond & A.S. Szalay, Astrophys. J. 274 (1983) 443. 
[53] P.J.E. Peebles, Science 224 (1984) 1385. 
[54] N. Kaiser, Astrophys. J. 273 (1983) L17. 

[55] S.D.M. White, M. Davis k C. Frenk, Mon. Not. R. Astr. Soc. 209 (1984) 27p. 

[56] Y. Zeng & S.D.M. White, Astrophys. J. 374 (1991) 1. 

[57] E. Braun, A. Dekel & P.R. Shapiro, Astrophys. J. 328 (1988) 34. 

[58] J. Centrella et al. , Astrophys. J. 333 (1988) 333. 

[59] R. Cen & J.P Ostriker, Astrophys. J. 399 (1992) 331. 

[60] P.J.E. Peebles, Astrophys. J. 263 (1982) LI. 

[61] J.R. Bond & G. Efstathiou, Astrophys. J. 285 (1984) L45. 

[62] M. Davis, G. Efstathiou, C.S. Frenk & S.D.M. White, Astrophys. J. 292 (1985) 371. 
[63] J.P. Ostriker, Ann. Rev. Astron. Astrophys. 31 (1993) 689. 

[64] G. Jungman, M. Kamionkowski & K. Griest, Phys. Rept. 267 (1996) 195 [arXiv:hep- 
ph/9506380]. 

[65] A.D. Linde, Particle Physics and Inflationary Cosmology (Harwood Academic, 1990). 

[66] A.R. Liddle h D.H. Lyth, Cosmological Inflation and Large-Scale Structure (Cambridge 
Univ Press, 2000). 

[67] R.K. Sachs & A.M. Wolfe, Astrophys. J. 147 (1967) 73. 

[68] R.H. Brandenberger, N. Kaiser, D.N. Schramm & N. Turok, Phys. Rev. Lett. 59 (1987) 
2371. 

16 



[69] E. Bertschinger k P.N. Watts, Astrophys. J. 328 (1988) 23. 

[70] J. Gratsias, R.J. Scherrer, G. Steigman k J.V. Villumsen, Astrophys. J. 405 (1993) 30. 

[71] P.J.E. Peebles, in The Origin and Evolution of Galaxies, ed. B. J. T. Jones (Reidel, 1983) 
p.143. 

[72] N. Sugiyama, M. Sasaki k K. Tomita, Astrophys. J. 338 (1989) L45. 

[73] E.L. Wright et al, Astrophys. J. 396 (1992) L13. 

[74] M. Davis, F.J. Summers k D. Schlegel, Nature 359 (1992) 393. 

[75] A.N. Taylor k M. Rowan-Robinson, Nature 359 (1992) 396. 

[76] D. Scott, J. Silk k M.J. White, Science 268 (1995) 829. 

[77] Q. Shan k F.W. Stecker, Phys. Rev. Lett. 53 (1984) 1292. 

[78] R.K. Schaefer, Q. Shan k F.W. Stecker, Astrophys. J. 347 (1989) 575. 

[79] A. Van Dalen k R.K. Schaefer, Astrophys. J. 398 (1992) 33. 

[80] A. Klypin, J. Holtzman, J. Primack k E. Regos, Astrophys. J. 416 (1993) 1. 

[81] Y.P. Jing, H.J. Mo, G. Borner k L.Z. Fang, Astron. Astrophys. 284 (1994) 703. 

[82] CP. Ma k E. Bertschinger, Astrophys. J. 429 (1994) 22. 

[83] D.Y. Pogosian k A. A. Starobinsky, Astrophys. J. 447 (1995) 465. 

[84] A.R. Liddle, D.H. Lyth, R.K. Schaefer, Q. Shah k P.T. Viana, Mon. Not. Roy. Astron. 
Soc. 281 (1996) 531. 

[85] Q. Shan k R.K. Schaefer, arXiv:hep-ph/9612478. 

[86] CP. Ma k E. Bertschinger, Astrophys. J. 434 (1994) L5. 

[87] J.R. Primack, J. Holtzman, A. Klypin k D.O. Caldwell, Phys. Rev. Lett. 74, 2160 (1995). 

[88] C. Athanassopoulos et al. [LSND Collaboration], Phys. Rev. Lett. 75 (1995) 2650. 

[89] Y. Fukuda et al. [Kamiokande Collaboration], Phys. Lett. B 335 (1994) 237. 

[90] M.J. White, D. Scott, J. Silk k M. Davis, Mon. Not. Roy. Astron. Soc. 276 (1995) L69. 

[91] G.G. Ross k S. Sarkar, Nucl. Phys. B 461 (1996) 597. 

[92] J.A. Adams, G.G. Ross k S. Sarkar, Phys. Lett. B 391 (1997) 271. 

[93] A.G. Riess et al. [Supernova Search Team Collaboration], Astron. J. 116 (1998) 1009. 

[94] S. Perlmutter et al. [Supernova Cosmology Project Collaboration], Astrophys. J. 517 (1999) 
565. 

[95] N.A. Bahcall, J. P. Ostriker, S. Perlmutter k P.J. Steinhardt, Science 284 (1999) 1481. 
[96] V. Sahni k A.A. Starobinsky, Int. J. Mod. Phys. D 9 (2000) 373. 
[97] O. Elgaroy et al, Phys. Rev. Lett. 89 (2002) 061301. 
[98] L.F. Abbott k M.B. Wise, Phys. Lett. B 135 (1984) 279. 

17 



[99] W. Hu, D. Scott, N. Sugiyama & M.J. White, Phys. Rev. D 52 (1995) 5498. 



[100 
[101 
[102 
[103 
[104 
[105 
[106 
[107 
[108 
[109 
[110 
[111 
[112 
[113 
[114 
[115 
[116 

[117; 
[us; 

[119 
[120 
[121 
[122 
[123 
[124 

[125 
[126 
[127 
[128 
[129 



E. Bertschinger, arXiv:astro-ph/9506070. 

U. Seljak & M. Zaldarriaga, Astrophys. J. 469 (1996) 437. 
D.H. Lyth, Phys. Rev. Lett. 78 (1997) 1861. 

W. Hu & S. Dodelson, Ann. Rev. Astron. Astrophys. 40 (2002) 171. 

W. Hu, N. Sugiyama & J. Silk, Nature 386 (1997) 37. 

P. de Bernardis et al [Boomerang Collaboration], Nature 404 (2000) 955. 

A. Balbi et al, Astrophys. J. 545 (2000) LI [Erratum-ibid. 558 (2001) L145]. 

S. Dodelson, E. Gates k A. Stebbins, Astrophys. J. 467 (1996) 10. 

X.M. Wang, M. Tegmark & M. Zaldarriaga, Phys. Rev. D 65 (2002) 123001. 

S. Hannestad, Phys. Rev. D 66 (2002) 125011. 

W. Hu, D.J. Eisenstein & M. Tegmark, Phys. Rev. Lett. 80 (1998) 5255. 
S. Buries, K.M. Nollett & M.S. Turner, Astrophys. J. 552 (2001) LI. 

C. L. Bennett et al, Astrophys. J. 583 (2003) 1. 

D. N. Spergel et al, arXiv:astro-ph/0302209. 

A. Pierce & H. Murayama, arXiv:hep-ph/0302131. 

C. Giunti, arXiv:hep-ph/0302173. 

G. Bhattacharyya, H. Pas, L. Song & T.J. Weiler, arXiv:hep-ph/0302191. 

H. V. Klapdor-Kleingrothaus, A. Dietz & I.V. Krivosheina, Found. Phys. 32 (2002) 1181. 
S. Dodelson et al [SDSS Collaboration], Astrophys. J. 572 (2001) 140. 

F. Hoyle & R.J. Tayler, Nature 203 (1964) 1108. 
P.J.E. Peebles, Phys. Rev. Lett. 16 (1966) 411. 

V.F. Shvartsman, Pisma Zh. Eksp. Teor. Fiz. 9 (1969) 315 [JETP Lett. 9 (1969) 184]. 

G. Steigman, D.N. Schramm & J.E. Gunn, Phys. Lett. B 66 (1977) 202. 

D. Denegri, B. Sadoulet & M. Spiro, Rev. Mod. Phys. 62 (1990) 1. 

J.M. Yang, M.S. Turner, G. Steigman, D.N. Schramm & K.A. Olive, Astrophys. J. 281 
(1984) 493. 

G. Steigman, K.A. Olive, D.N. Schramm & M.S. Turner, Phys. Lett. B 176 (1986) 33. 

J.R. Ellis, K. Enqvist, D.V. Nanopoulos & S. Sarkar, Phys. Lett. B 167 (1986) 457. 

G. Steigman, K.A. Olive k D.N. Schramm, Phys. Rev. Lett. 43 (1979) 239. 

V. Barger, P. Langacker & H.S. Lee, arXiv:hep-ph/0302066. 

J. Bernstein, L.S. Brown & G. Feinberg, Rev. Mod. Phys. 61 (1989) 25. 

18 



[130] S. Sarkar, Rept. Prog. Phys. 59 (1996) 1493. 

[131] K.A. Olive, D.N. Schramm, G. Steigman k T.P. Walker, Phys. Lett. B 236 (1990) 454. 

[132] T.P. Walker, G. Steigman, D.N. Schramm, K.A. Olive k H.S. Kang, Astrophys. J. 376 
(1991) 51. 

[133] C.J. Copi, D.N. Schramm k M.S. Turner, Science 267 (1995) 192. 

[134] N. Hata, R.J. Scherrer, G. Steigman, D. Thomas, T.P. Walker, S. Bludman k P. Lan- 
gacker, Phys. Rev. Lett. 75 (1995) 3977. 

[135] P.J. Kernan k S. Sarkar, Phys. Rev. D 54 (1996) 3681. 

[136] S. Sarkar, arXiv:astro-ph/0205116. 

[137] B.D. Fields k K.A. Olive, Astrophys. J. 506 (1998) 177. 
[138] K.A. Olive k E. Skillman New Astron. 6 (2001) 119. 
[139] D. Sauer k K. Jedamzik, arXiv:astro-ph/0104392. 
[140] Y.I. Izotov et al. , Astrophys. J. 527 (1999) 757. 

[141] M. Peimbert, A. Peimbert k M.T. Ruiz, Astrophys. J. 541 (2000) 688. 
[142] A. Peimbert, M. Peimbert k V. Luridiana, Astrophys. J. 565 (2002) 668. 
[143] S.G. Ryan, J.E. Norris k T.C. Beers, Astrophys. J. 523 (1999) 654. 
[144] E. Vangioni-Flam et al. , New Astron. 4 (1999) 245. 

[145] S.G. Ryan, T.C. Beers, K.A. Olive, B.D. Fields k J.E. Norris, Astrophys. J. 530 (2000) 
L57. 

[146] S. Vauclair k C. Charbonnel, Astrophys. J. 502 (1998) 372. 

[147] M.H. Pinsonneault, T.P. Walker, G. Steigman k V.K. Narayanan, Astrophys. J. 527 
(1999) 180. 

[148] J. Linsky, Space Sci. Rev. 84, 285 (1998). 

[149] A. Songaila, L.L. Cowie, C. Hogan k M. Rugers, Nature 368 (1994) 599 . 
[150] J.K. Webb et al. , Nature 388 (1997) 250. 

[151] S. Buries, D. Kirkman k D. Tytler, Astrophys. J. 519 (1999) 18. 
[152] D. Kirkman et al, Astrophys. J. 559 (2001) 23. 
[153] J.M. O'Meara et al. , Astrophys. J. 552 (2001) 718. 

[154] D. Kirkman, D. Tytler, N. Suzuki, J.M. O'Meara k D. Lubin, arXiv:astro-ph/0302006. 
[155] S. D'Odorico, M. Dessauges-Zavadsky k P. Molaro, Astron. Astrophys. 368 (2001) L21. 
[156] M. Pettini k D. Bowen, Astrophys. J. 560 (2001) 41. 

[157] S.A. Levshakov, M. Dessauges-Zavadsky, S. D'Odorico k P. Molaro, Astrophys. J. 565 
(2002) 696. 

19 



[158 
[159 

[160 
[161 
[162 
[163 
[164 
[165 

[166 
[167 
[168 
[169 
[170 
[171 
[172 
[173 
[174 
[175 
[176 

[177; 
[178; 

[179 

[i8o; 

[181 
[182 

[183; 

[184 
[185 



G. Sonneborn et al. , Astrophys. J. 545 (2000) 277. 

D. Kirkman, D. Tytler, S. Buries, D. Lubin k J.M. O'Meara, Astrophys. J. 529(2000) 
655. 

T.M. Bania, R.T. Rood k D.S. Balser, Nature 415 (2002) 54. 

D. Galli, L. Stanghellini, M. Tosi k F. Palla, Astrophys. J. 477 (1997) 218. 
C. Charbonnel, Space Sci. Rev. 84(1998) 199. 

R.H. Cyburt, B.D. Fields k K. A. Olive, Astropart. Phys. 17 (2002) 87. 

J.A. Adams, G.G. Ross k S. Sarkar, Nucl. Phys. B 503 (1997) 405. 

J. Barriga, E. Gaztanaga, M.G. Santos k S. Sarkar, Mon. Not. Roy. Astron. Soc. 324 
(2001) 977. 

S.L. Bridle, A.M. Lewis, J. Weller k G. Efstathiou, arXiv:astro-ph/0302306. 

E. Lisi, S. Sarkar k F.L. Villante, Phys. Rev. D 59 (1999) 123520. 

S. Buries, K.M. Nollett, J.N. Truran k M.S. Turner, Phys. Rev. Lett. 82 (1999) 4176. 
P. Di Bari, Phys. Rev. D 65 (2002) 043509. 
K. N. Abazajian, arXiv:astro-ph/0205238. 
S. Hannestad, arXiv:astro-ph/0302340. 

R.H. Cyburt, B.D. Fields k K.A. Olive, arXiv:astro-ph/0302431. 
E. Pierpaoli, arXiv:astro-ph/0302465. 

P. Crotty, J. Lesgourgues k S. Pastor, arXiv:astro-ph/0302337. 

H. -S. Kang k G. Steigman, Nucl. Phys. B372, 494 (1992). 

C. Lunardini k A.Y. Smirnov, Phys. Rev. D 64 (2001) 073006. 

A.D. Dolgov, S.H. Hansen, S. Pastor, S.T. Petcov, G.G. Raffelt k D.V. Semikoz, Nucl. 
Phys. B 632 (2002) 363. 

Y.Y. Wong, Phys. Rev. D 66 (2002) 025015. 

K.N. Abazajian, J.F. Beacom k N.F. Bell, Phys. Rev. D 66 (2002) 013008. 

J.H. Applegate, C.J. Hogan k R.J. Scherrer, Phys. Rev. D 35 (1987) 1151. 

K. Kainulainen, H. Kurki-Suonio k E. Sihvola, Phys. Rev. D59, 083505 (1999); 
K. Jedamzik k J.B. Rehm, Phys. Rev. D64, 023510 (2001). 

J.B. Rehm k K. Jedamzik, Phys. Rev. Lett. 81, 3307 (1998); 
H. Kurki-Suonio k E. Sihvola, Phys. Rev. D62, 103508 (200). 

S. Dimopoulos, R. Esmailzadeh, G.D. Starkman k L.J. Hall, Astrophys. J. 330 (1988) 
545; 

K. Jedamzik, Phys. Rev. Lett. 84 (2000) 3248. 

N. Arkani-Hamed, S. Dimopoulos, G.R. Dvali k J. March-Russell, Phys. Rev. D 65 (2002) 
024032. 

K. Abazajian, G.M. Fuller k M. Patel, arXiv:hep-ph/0011048. 



20 



